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Abstract 

We study universality and scaling properties of RG gauge actions (Iwasaki 
and DBW2). In the first part we consider the critical temperature T c and com- 
pute the reference energy scale r for critical couplings /3 C corresponding to N t = 
3, 4, 6, 8. The universality of T c ro between Iwasaki and Wilson action is confirmed 
and the scaling behavior of the Iwasaki action is found to be better than the 
one for the Wilson action. The results for the DBW2 action show larger lattice 
artefacts. A continuum value T c r = 0.7498(50) is extracted. We compute also 
the glueball masses for the states ++ and 2 ++ , investigate the scaling of m ++r Q 
and 777,2++ ro and point out practical problems which are due to the violation of 
positivity present in the RG actions. 



1 Introduction 



In view of the next unquenched lattice simulations, big efforts were devoted in the 
last years to improve both fermionic and gauge actions. 

Although the lattice artefacts for the standard gauge action start at 0(a 2 ) and 
hence one expects that these are somehow less relevant than the ones induced by 
the fermionic part, it turns out that the gauge term plays an important role, and 
the question which is the most convenient gauge action to adopt has been often 
addressed. 

The purpose of adopting alternative actions is not only to improve the scaling 
behavior; also features related to chiral symmetry were investigated. In particular, 
RG actions (Iwasaki, DBW2) have been suggested as good candidates to be used 
in the next simulations on Ginsparg- Wilson/domain wall fermions; interesting 
characteristics were observed, such as the suppression of small instantons and 
dislocations and a possible remedy of the problem of residual chiral symmetry 
breaking for domain wall fermions [1,2]. Other authors [3], however, pointed out 
possible problems connected to these actions and proposed an alternative method 
to reduce residual-mass effects that works very well also with the Wilson plaquette 
action. 

This increasing interest and discussion in improved gauge actions motivates more 
investigations into their properties, starting from the basic ones, like universality 
and scaling behavior. 

There are in principle several quantities that one can use to quantify the lattice 
artefacts and to test universality by comparing the results with the plaquette 
action known in the literature; in particular in this work we will consider the 
critical deconfining temperature T c and the glueball masses for the states ++ and 
2++. 

In a previous work [4], the scaling of the ratio T c /^/a has been investigated and 
lead for the Iwasaki action to a continuum limit which is in disagreement with 
the one obtained with the Wilson action. More likely than a possible violation of 
universality, this fact is related to the definition of the string tension a. 
We will perform the same analysis but using the scale r ~ 0.5fm [5] instead of a. 
It is well known that the extraction of the string tension from the static potential 
at finite quark separations is problematic and leads to systematic errors, which 
are difficult to control. It is preferable to define the scale by using intermediate 
distance properties. 

After computing r^/a at several lattice spacing for Iwasaki and DBW2 actions we 
will perform a continuum extrapolation of the quantity T c r , discuss universality 
and compare the scaling violations for the different actions. 
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We will point out possible problems that can occur by adding irrelevant operators 

to the plaquette action, as the violation of physical positivity. 

The second part of the work is dedicated to the evaluation of the glueball masses. 

In particular the ++ mass is a promising observable to discuss the lattice artefacts, 

since for the Wilson action large scaling violations have been observed. 

Despite the large errors and the difficulties in the calculations, we will discuss our 

results, draw indicative conclusions and remark on open questions. 

2 Improved gauge actions 

In this work we will consider SU(3) Yang-Mills lattice actions formulated on the 
basis of renormalization group (RG) considerations. The most popular exam- 
ples are the Iwasaki [6] and DBW2 [7,8] actions, which are restricted to a two- 
parameter space and include planar rectangular (1 x 2) loops in addition to the 
usual plaquette term, 

S = ^E ( c °E i 1 - ! Re ^ x V)} + c i E i 1 - ! ReW ^ x2 (*)} ) , (2.i) 

X \ ^1<V fj,,l/ / 

with the normalization condition Co = 1 — 8c\. 

The coefficient C\ in eq. (2.1) takes different values for various choices of alternative 
actions 

{—1/12 Symanzik, tree level impr. 
-0.331 Iwasaki, RG (2.2) 
-1.4088 DBW2, RG 

where we included also the 0(a 2 ) tree-level improved Symanzik action [9], [10]. 
First of all one can notice that the strength of the rectangular loops for RG 
improved actions is significantly larger than what is needed in order to cancel the 
0(a 2 ) effects at tree level. We will later discuss this further. 
We will also consider the gauge FP action [11] for comparison, in particular the 
study on the lattice artefacts performed in [12]. 

3 The critical temperature T c 

It is well known [13,14] that pure Yang-Mills theory undergoes a first order phase 
transition at some finite temperature T c . On the lattice, the critical temperature 
is determined by evaluating the critical coupling f3 c 

£ = N t a(p c ), (3.1) 
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N t 


(3 C Wilson [16] 


Iwasaki [4] 


DBW2 [17] 


3 




2.1551(12) 


0.75696(98) 


4 


5.6925(2) 


2.2879(11) 


0.82430(95) 


6 


5.8941(5) 


2.5206(24) 


0.9636(25) 


8 


6.0624(12) 


2.7124(34) 




12 


6.3380(23) 







Table 1: The critical coupling for Wilson, Iwasaki and DBW2 actions. 

where N t is the number of lattice points in the time-like direction with periodic 
boundary conditions. The extension in the space-like directions is supposed to be 
infinite. 

There are several methods for determining f3 c ; for example, one can locate the 
peak in the Polyakov loop susceptibility. 

In addition to its intrinsic importance as a fundamental non-perturbative predic- 
tion, T c provides also a useful quantity to study the lattice artefacts for different 
gauge actions and to test universality. 

For known values of (3 C at given N t for different actions one can refer to [15,16,4,17, 
12]; for Wilson, Iwasaki and DBW2 actions, these values are collected in Table 1. 



action 


N t 




Wilson [16] 


oo 


0.630(5) 


Iwasaki [4] 


oo 


0.651(12) 


DBW2 [17] 


oo 


0.627(12) 


Sym. tree level [16] 


oo 


0.634(8) 


1 loop tadpole impr. [18] 


oo 


0.659(8) 


FP [12] 


4 


0.624(7) 



Table 2: Results for the deconfming temperature in units of the string tension from different 
actions. The continuum extrapolations are taken from [19]. 

The available results in the literature are mostly expressed in terms of the 
string tension a; the results for the quantity T c /^fa are reported in Table 2, 
where the continuum extrapolations are from [19]. Figure 1 collects the results 
for Wilson, Iwasaki, DBW2, Symanzik tree level and FP action, and was taken 
from [12], where the latest evaluation of T c /y/a (FP action) has been performed. 

The first observation from Fig. 1 is that for this specific quantity the dis- 
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Figure 1: T c /y/a as function of for different actions, from [12]. By "RG improved" is 

meant here the Iwasaki action. 

cretization effects appear not to be very significant and hence is not possible to 
arrive at precise conclusions about the lattice artefacts for different actions. 

Furthermore, one notices that for the Iwasaki and the Wilson actions a dif- 
ference of order 2a in the continuum results is observed. The most drastic expla- 
nation for this discrepancy could be a violation of universality but this scenario 
seems unrealistic; a more natural explanation is that the string tension is difficult 
to determine and systematic errors due to this were not included in all calcula- 
tions. It is preferable to use r [5] to reliably set the scale. 

Let us recall that the string tension is related to the properties of the force be- 
tween static quarks for distances r — > oo, and hence its evaluation at finite r can 
in principle contain large systematic errors. 

On the contrary, r ~ 0.5fm is extracted from the force at intermediate distances 
and can be evaluated very precisely. For the Wilson action this quantity has been 
evaluated in the coupling range 5.7 < (3 < 6.92 [20,21]. 



4 Evaluation of ro/a for RG actions 

For the Wilson action, the values of ro/a corresponding to the critical couplings 
at different N t can be easily obtained by the parametrization formula in [20]. 
For the Iwasaki and DBW2 actions there was up to now no precise evaluation of 
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ro/ a and we performed new numerical simulations with this purpose. 
In our computation we followed essentially the procedure adopted in [5,20] for the 
plaquette action. We applied the smearing procedure [22] to the spatial links and 
the multi-hit method [23] to the time-like links for the variance reduction. 
For each spatial separation r we constructed Wilson loop correlation matrices 

c lrn {t) = (tt [v t (o, A)v(A, A + td)vX(tb, A + m)v\o, to) } ^ (4.1) 

with l,m — 1, ...,M; Vi(x, y) indicates the product of link variables connecting x 
and y at smearing level / in a spatial direction, while V(x, y) is the product of 
time-like links between x and y after the application of the multihit procedure. 
In particular we have chosen M — 4, and the number of smearing iterations 
for each level has been determined with the same criterion as for the Wilson 
action [20]. 

We adopted a hybrid algorithm with N or over-relaxation steps per heat-bath step 
[24,25] and increased N or with (3 according to 

N or « 1.5(r /a). (4.2) 

The simulation parameters are reported in Table 3. For the DBW2 action, besides 
the three values of f3 c known in the literature (Table 1), we decided to evaluate 
r /a also for a larger f3 = 1.04, which should roughly correspond to (5 — 6 for the 
Wilson action and has been used in quenched simulations [2]. 

We started our analysis by observing the time-dependence of the effective 
potential evaluated from the diagonal elements of the correlation matrices. 
Figure 2 shows the effective potential 

^(r) = -ln(^)) (4.3) 

evaluated for I — m — 2, which corresponds to what was estimated to be the 
optimal smearing for the Wilson action in [26]. 

It has been pointed out in [27] (see app. A) that actions containing terms 
in addition to the plaquette violate the physical positivity, and hence one expects 
also negative contributions to the spectral decomposition of correlation functions. 
These negative contributions for the RG actions are evident in Fig. 2 and, as 
expected from our discussion in app. A, they are invisible starting from a certain 
t which depends on the coefficient c\ of eq. (2.2). 

The most important consequence of this phenomenon regards the applicability of 
the variational method [28,29], which is mathematically founded on the positivity 
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Figure 2: The effective potential as function of t for r w r (r /a ~ 4 — 5) for the Wilson 
action and for the RG-improved Iwasaki and DBW2 action. 



of the correlation matrix Ci m (t) at a certain small t — to; this condition is verified 
only for t ^> t min , but on the other hand t can not be arbitrarily large because 
the statistical errors increase exponentially with t and make the inversion of Ci m 
impracticable. 

One observes however a quite satisfactory plateau in the effective potential, start- 
ing at sufficiently large t. In that sense the violation of physical positivity does 
not represent a real trouble for these correlation functions, but can become quite 
problematic for example for the extraction of the glueball masses (see sec. 7). 
We decided to extract the potential from eq. (4.3) at t/a = (3 — 4) without ap- 
plying the variational method. The systematic error was estimated by taking the 
difference between this value and what one would obtain extracting the potential 
at (t + a). We linearly added systematic and statistical errors; at small r the 
total uncertainty is dominated by the systematic one, while at large distances the 
situation is in general reversed. 



4.1 The force at tree level 

A tree-level study of the force for the different actions can furnish important hints 
on how the continuum limit is approached, in particular in the region of small 
couplings. 
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Figure 3: The force at tree level for several actions. 



Explicitly, for the action eq. (2.1) the tree-level force is given by [9] 

F(r') = y(r) '^ (r ~ a) = (4.4) 

_4£o / d A k cos(rfci/a)-cos((r-a)fci/a) . r\( A 2\ _ 

3 a J ^ (2^)3 4(E J 3 =lS m 2 (fc i /2)-4 Cl E J 3 . 1 sm 4 (fc J /2)) + U ^0« ) 

= F tree (r') + 0( g y), 

where — 6/(3. 

In Fig. 3 the quantity ^--F tree (r n ), with r n = r — a/2 is plotted as function of 
(a/r n ) 2 for different actions. One can notice that the RG actions, in particular 
the DBW2, show at tree level large lattice artefacts at small (a/r'). In that sense, 
at tree level the RG actions are "over-corrected" and introduce lattice artefacts 
of the same order or even larger than what is expected with the usual plaquette 
action. This fact indicates that the continuum extrapolation should be considered 
with great care, because unless {a/r') is very small, one can not be sure to be in 
the region where the leading discretization errors are quadratic in a. 
Beside the naive definition of the force 

F(r') = [V(r) - V(r - a)]/ a, (4.5) 

with r' = r n , we followed the procedure reported in [5] by introducing a tree- level 
improved definition with r' = n such that 

Ftreein) = | ^ , (4.6) 
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with no lattice artefacts at tree level. For the Wilson action this choice turned 
out to improve substantially the scaling behavior of the force [5,20]. 
From eq. (4.4) one obtains 

(4 7 rr I 2 )- 1 = -[G(r,0,0) - G(r - a,0,0)]/a, (4.7) 

where 

r(iA — if J^t ll|=iCos(r^/a) x 

a J ^ (2^)3 4(E'U sm 2 (fc J /2)-4c 1 £* =1 sin 4 (fc i /2)) ' 

is the scalar free lattice propagator associated to the action eq. (2.1). We computed 
ri by solving the integral eq. (4.8) numerically. 

Iwasaki action: 



L/a 


P 


ni 


N 

iy or 


N 

1 v meas 


8 


2.1551 


0,2,4,6 


3 


20000 


12 


2.2879 


0,4,9,13 


4 


4000 


24 


2.5206 


0,12,25,37 


8 


645 


32 


2.7124 


0,18,36,54 


9 


370 


DBW2 action : 


L/a 


P 


ni 


N or 


Nmeas 


10 


0.75696 


0,2,4,6 


3 


12000 


12 


0.8243 


0,4,9,13 


4 


6000 


16 


0.9636 


0,10,20,30 


8 


800 


24 


1.04 


0,18,36,54 


9 


220 



Table 3: Simulation parameters for the Iwasaki and DBW2 actions for the evaluation of ro/a. 
L is the lattice extension (the spatial and temporal extensions are equal), ni represents the 
number of smearing iterations for each level. 

4.2 Results 

Our numerical results for the potential and the force at finite lattice spacing are 
collected in the tables 11 and 12 in App. B. 

Once the force has been evaluated, we extracted the value of r /a by using a local 
interpolation formula; the results are reported in Table 4. We adopted both the 
naive and the tree-level improved definition of the force. 

The first error contains the statistical uncertainty summed to the systematic 
one due to the interpolation of the force. The second error is the systematic 
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Iwasaki action : 



p 


ro/a (r n ) 


To/ a (n) 


r /a (eq. (4.9)) 


2.1551 


2.311(5)(9) 


2.320(6)(9) 


2.333 


2.2879 


3.026(4)(3) 


3.026(5)(1) 


3.021 


2.5206 


4.535(6)(4) 


4.511(8)(1) 


4.514 


2.7124 


6.020(15)(25) 


5.999(15)(19) 


5.983 



DBW2 action : 



p 


ro/a (r n ) 


r /a (n) 


0.75696 


2.430(5)(20) 


2.225(4)(11) 


0.8243 


3.129(23)(1) 


3.036(17)(4) 


0.9636 


4.606(13)(17) 


4.556(17)(20) 


1.04 


5.500(29)(7) 


5.452(26)(8) 



Table 4: Results for ro/a evaluated at different (i = (3 C for Iwasaki and DBW2 actions, using 
the naive definition of the force or the tree-level improved eq. (4.7). For the Iwasaki action, the 
fourth column is obtained from the fits eq. (4.9). 

uncertainty coming from different choices of t in the effective potential eq. (4.3). 
One can notice that for the DB W2 action the choice of r n or r\ in the definition 
of the force leads to results which can be quite different from each other, above 
all at small r /a; we expected this feature by investigating the force at tree level. 
This ambiguity will make the discussion of the lattice artefacts difficult, because 
the possible conclusions will depend on which definition of the force one has used 
and not on intrinsic properties of the action. 

For the Iwasaki action the results obtained through the two definitions are not 
significantly different. 

4.3 Parametrization of r /a 

Following the strategy of [20] , one can attempt a phenomenological parametriza- 
tion of r /a in the range of couplings under consideration. 

For the Iwasaki action the four values of r (obtained by adopting r{) were fitted 
in the form 

ln(a/r ) = d + c 2 (/3 - 3) + c 3 (/3 - 3) 2 , (4.9) 
yielding the numerical results 

ci = -2.1281, c 2 = -1.0056, c 3 = 0.6041. (4.10) 
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Figure 4: Paramctrization of r /a for the Iwasaki action, eq. (4.9), using the tree-level improved 
definition of the force. 

in the range 2.1551 < (3 < 2.7124. 

The results and the fit formula are shown in Fig. 4; the accuracy is about 0.6% 
at (3 = 2.1551 and 0.8% at (3 = 2.7124. 

For the DBW2 action we had to use a four-parameter representation 

ln(a/r ) = d ± + d 2 {(3 - 1) + d z {(3 - l) 2 + d 4 (f3 - l) 3 , (4.11) 

with 

= -1.6007, d 2 = -2.3179, d 3 = -0.8020, d A = -19.8509, (4.12) 

for the range 0.75696 < (3 < 1.04, where the results always refer to the tree-level 
improved definition of the force. The parametrization is plotted in Fig. 5. 
Note that in [2] the value r /a = 5.24(3) for the DBW2 action at (3 = 1.04 is 
quoted. Our value differs about 3% from that one. 

5 Scaling of T c tq 

Once r /a at the given couplings is known, one can finally consider the renormal- 
ized quantity 

r oT c = (5.1) 
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P 

Figure 5: Parametrization of r /a for the DBW2 action, eq. (4.11), using the tree-level im- 
proved definition of the force. 

The results for Iwasaki and DBW2 actions are given in Table 5, together with 
the values obtained using the Wilson action. Also in this case, we show both 
the results obtained with the naive and with the tree-level improved definition of 
the force. The error in T c tq is the quadratic sum of the error for r$/a and the 
uncertainty in f3 c , which can be translated into an uncertainty in r by using the 
parametrization formulas eq. (4.9), eq. (4.11). In our evaluations the error for f3 c 
and the uncertainty in r are roughly of the same order. 

We expect that the leading lattice artefacts are of order a 2 , such that the contin- 
uum limit is approached in the following way 

T c r = T c r \ a=0 + s ■ (aT c ) 2 + 0(aT c ) 4 . (5.2) 

The results for T c tq, together with the continuum extrapolation, are shown in 
Fig. 6. 

For the Iwasaki action there is no appreciable difference between the results 
obtained with r n and rj and in both cases the data show better scaling properties 
in comparison to the Wilson action. Furthermore, the value obtained at N t — 8 is 
in full agreement with the continuum result evaluated through the Wilson action 
and hence the universality is confirmed; this supports the conclusion that the 
disagreement observed in T C J \fa is indeed due to the difficulty in evaluating the 
string tension, particularly at small lattice spacings, and it is necessary to set the 
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Figure 6: T c r for different actions. The x coordinates were slightly shifted for clarity. 
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Figure 7: Continuum extrapolation of T c r for the Iwasaki and Wilson action, using the 
constrained fit eq. (5.2). 
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c u 

N t 


Wilson 


Iwas.frr,) 


Iwas.(ri) 


DBW2(r n ) 


DBW2(ri) 


3 




0.7703(50) 


0.7733(53) 


0.8100(90) 


0.7417(90) 


4 


0.7194(22) 


0.7565(23) 


0.7565(23) 


0.7822(65) 


0.7590(56) 


6 


0.7388(30) 


0.7558(33) 


0.7518(31) 


0.7676(68) 


0.7593(76) 


8 


0.7445(33) 


0.7525(60) 


0.7499(55) 






12 


0.7464(44) 











Table 5: Results for T c r . 
scale through a more reliable quantity. 

Also for the DBW2 action the scaling properties are improved, although only by 
adopting 77 instead of r n , so that it is more difficult to make a statement about 
the lattice artefacts in this case. 

A constrained fit of the form eq. (5.2) including the points with N t > 6 for Iwasaki 
and Wilson actions yields the continuum result (Fig. 7) 

T c r = 0.7498(50). (5.3) 

At N t = 6 the Wilson action shows scaling violations for r T c of about 1.5%, while 
they are 0.3% for the Iwasaki action. 

For N t — 4 the discretization errors for the Wilson action increase to 4%, while 
for the Iwasaki action they remain very small (0.6%). 

In Fig. 6 we included also the results obtained with the FP action [12], which also 
show a good scaling within 1% even on coarse lattices corresponding to N t = 3, 2. 
One has however to mention that for those lattices the determination of r /a 
contains large systematic uncertainties, as pointed out by the authors. 

6 Scaling of a q q(//) 

Another interesting observable that can be used to test scaling violations is the 
dimensionless coupling a qq -(/j) obtained from the force. For this purpose we com- 
pared our present determination of the force at finite lattice spacing with the 
results of [21], where the continuum extrapolation has been performed in the re- 
gion 0.05fm < r < 0.8fm. 

We point out that we only determined the on-axis potential and hence we can not 
investigate violations of rotational invariance which would require the evaluation 
of off- axis quantities. 
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Figure 8: a q q at finite lattice spacing for Iwasaki and DBW2 action compared with the con- 
tinuum result. The solid line represents the 3- loop RG perturbative prediction of the running 
coupling; the dashed lines correspond to its uncertainty [30] . 



The coupling a qq - is denned in terms of the force by the simple relation 

a m (pi) = f F(r)r 2 , n=\. (6.1) 

Figure 8 shows a; qq - in the continuum limit and the results obtained with the 
Iwasaki and DBW2 actions at the largest (3 at our disposal. For the Iwasaki action 
no appreciable difference in the results obtained with r n and r\ can be seen, while 
for the DBW2 action the discrepancy becomes large at small distances. 
At large enough distances one obtains a good scaling in the coupling, and one does 
not observe scaling violation within the statistical errors. At small r/r one sees 
deviations from the continuum limit, as one can observe in Fig. 9, where only the 
short distance region is considered. At r ~ OAtq and a ~ 0.09fm, the deviations 
can be estimated to about 2% for the Iwasaki action; for the DBW2 action they 
amount to 4% if one uses r n to define the force, and reach even 40% by employing 
r\. This fact shows that the adoption of a tree level improved definition of the 
force does not guarantee success in reducing the lattice artefacts; in particular for 
RG actions, which are over-corrected at tree level, one should always check the 
scaling violations for different observables in order to obtain safe statements. 

Note that for the DBW2 action the lattice artefacts in T c r are large when 
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Figure 9: a q q at finite lattice spacing for Iwasaki and DBW2 action compared with the con- 
tinuum results in the short distance region. 

one uses r n to define the force. Thus there is not one definition of the force at 
finite a, which has reasonably small lattice artefacts for both T c r and a q q. 

7 Glueball masses 

For the computation of the glueball masses we followed the method proposed 
by [12]. 

We decided to concentrate on the states J PC = ++ , 2 ++ by measuring the masses 
in the irreducible representations Af + , E ++ and T 2 ++ of the cubic group. A\ 
corresponds to the J = state; in the continuum limit one expects that the 0(3) 
symmetry is restored and hence the doublet E and the triplet T 2 are degenerate 
and form together the J = 2 quintuplet. 

The m ++ is particularly interesting for our purpose to investigate the lattice 
artefacts on RG actions since these turn out to be quite sizeable for the Wilson 
action. 

We performed a test simulation using the standard Wilson action for (5 = 6.0 on 
a 16 4 lattice; here we measured all 22 loop shapes up to length 8 that can be 
build on the lattice (see [31] for the classification of the loops) and formed the 
wave functions corresponding to the representations R = Af + ,E ++ ,T^~ + . Our 
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timeslice observable at time t is 

dn 

^) = ^EE< fiRe ^^)' n = l,...,22, (7.1) 

where the coefficient c 1 ^ are taken from the literature [31,32], the sum over % 
indicates the sum over all d n orientations of a given shape n and K is a suitable 
normalization constant. 
Then we build the correlation matrices 1 

cm = {s*{t)sm)c = {s*{t)s*m - (sm(sm)- (7.2) 

The indices k, I assume (22 x M) values, where M is the number of smearing levels. 
We adopted the same smearing procedure [22] that we used for the computation 
of the Wilson loops. 

Starting from this large basis, we analyzed the signal/noise ratio of the dif- 
ferent operators in order to eliminate those which introduce large noise and to 
reduce the correlation matrix to a set of well measurable operators. 
We found general agreement in the classification of the "bad" operators with [12]. 
Our final choice for the operators to measure was then 

A+ + : #2, #5, #7, #8, #10, #12, #14 (7.3) 
E ++ : #2, #5, #7, #10, #12, #14 (7.4) 
T++ : #7, #12, #14, (7.5) 

where for E ++ and T 2 I " + channels we took respectively 2 and 3 orthogonal projec- 
tions for each shape. For each operator we then considered all M smearing levels. 
We assumed that this choice is reasonable for each value of the lattice spacing and 
also for the Iwasaki and DBW2 actions and we stress that we did not perform a 
systematic study on the features of the several operators. 

The simulation parameters for our new Monte Carlo simulations are listed in 
Table 6. Measurements were taken after a number of sweeps between 3 and 5. 

7.1 Analysis details 

As already explained in sec. 4, due to the violation of physical positivity for the 
RG actions (see app.A), the variational method [28,29] is mathematically not well 
founded, at least not at small time separations, where one would like to apply 

1 Notice that the vacuum subtraction is required only in the channel, since it has the 
same quantum numbers as the vacuum. 



16 



Iwasaki action 



p 


L 


ni 




2.2423 


10 


2,4,6,8 


12000 


2.2879 


12 


2,4,6,8 


16000 


2.5206 


16 


3,6,9,12 


8000 


DBW2 action 


P 


L 


ni 


^meas 


0.8342 


12 


2,4,6,8 


8000 


0.9636 


16 


3,6,9,12 


2500 



Table 6: Simulation parameters for the evaluation of the glueball masses for the Iwasaki and 
DBW2 actions. 

it. The statistical errors are indeed drastically increasing with t and already at 

t = 4a the signal in the correlation function eq. (7.2) is lost in noise. 

We decided to follow [12] by first solving the generalized eigenvalue problem 

C{h)v a {t u t„) = A a (ti, f )C(foK(fi, h), (7.6) 

with t = 0, t\ = a. Then we projected the correlation matrices to the space of 
eigenvectors corresponding to the N eigenvalues which satisfy the condition 2 

\p>e, P = l,...,N, (7.7) 

where e is an adjustable (small) parameter. So the reduced matrix is obtained by 

C%(t) = (vifaMCWv^to)), i,j = 1,...,N, (7.8) 

where the index R for the representation in now omitted. 

By choosing e appropriately in eq. (7.7) one hopes to get rid of most of the 
unphysical modes caused by negative and very small eigenvalues. 
Then one can apply the variational method to the reduced matrix 

C N (t)wp(t,to) = \p(t,t )C N (t )wp(t,t ). (7.9) 

The effective glueball mass can be read off directly from the largest eigenvalue 
corresponding to the lowest energy 

rneftf) = ~ (^Sy) ■ (7.10) 

2 In [12] the direct eigenvalues and eigenvectors of C(io) are considered, instead of the gener- 
alized one. We have tried both possibilities and found consistent results. 
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or, alternatively one can project again C to the subspace corresponding to the 
largest eigenvalue 



We tested that these two different evaluation yield results which are compatible 
within the statistical errors. 

We chose e in eq. (7.7) so that the reduced matrix had dimension between 2 and 
5. 

We applied the same procedure also for t — a, ti — 2a in eq. (7.6), but due to the 
fact that the statistical fluctuations are already quite large one has to start from 
the beginning from a reduced number of operators in order to be able to solve the 
generalized eigenvalue problem. 

We decided to extract the glueball masses by taking the effective mass at t = 3a 
for the small (3 regime and t = 4a for the large value of the coupling at our 
disposal. For the RG actions we performed numerical simulations up to L = 16, 
with a minimum lattice spacing a ~ O.lfm, which is quite large and hence we will 
not be able to perform a continuum extrapolation of the results. On the other 
hand, this is the regime of lattice spacings that have been used for simulations 
with dynamical fermions and for this reason it is desirable to obtain informations 
about the discretization errors in this region. 

The figures 10,11 show the effective masses in the A± + channel computed by 
using eq. (7.10) (filled squares) and eq. (7.12) (empty squares) for the Iwasaki 
and the DBW2 action, with to — 0. For the largest value of f3, we show also 
the results obtained with to — a (filled triangles). As for the potential, one 
notices also in this case the presence of negative contributions in the correlation 
functions due to unphysical states; the situation here is somehow more drastic 
since one has to discard for this reason the small t region but on the other hand 
the errors increase very rapidly and one is forced to extract the glueball mass at 
t = (3 — 4)a. The results can then be affected by systematic errors that can not be 
easily estimated. In the plots we have indicated with the dotted lines the range in 
which we decided to take the mass. The statistical errors were evaluated by using 
a jackknife procedure. 

Concerning the determination of m 2 ++ we observed that the signal for the E ++ 
channel is usually worse than for the T 2 ++ and the errors on the effective masses 
are very large already at t — 3a. For this reason we decided to use m T ++ as 
estimate of m 2 ++ at finite lattice spacing. 




(7.11) 




(7.12) 
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Figure 12 shows the effective masses for the T 2 ++ channel. The determination of 
the masses is more problematic that the Af + case; in particular for the largest 
coupling one can observe that at t — 4a the errors are too large to have a significant 
measurement of the mass. We decided hence to extract the mass at t — 3a, taking 
care that this value is compatible with the one at t — 4a within the statistical 
errors. 

7.2 Results 

In Table 8 the results for the masses in lattice units are reported. For the smallest 
P for the Iwasaki action, it was not possible to obtain a reliable evaluation of the 
mass in the T 2 ++ channel. 

The renormalized quantities r m G are reported in tables 9 and 10; due to the 
fact that the errors on am G are quite large (between 8% and 9% for Af + and 
between 10% and 20% for T^~ + ), the difference which arises by choosing 77 or r n 
is in this case much smaller than the total uncertainty; we have reported however 
both results in the table. For f3 = 2.2423 we had no direct measurement of r^/a 
and we made use of the interpolation formula eq. (4.9). 

Figure 13 shows the results for r m ++ as function of (a/r ) 2 , displaying only the 
results obtained with r n . For the comparison we included the results for FP ac- 
tion [12] and several calculations performed with the Wilson action [33,34,35]. 
The continuum values available in the literature are listed in Table 7 and have 
been taken from [12] for the FP action and from [36] for the Wilson action, where 
the results of [33,34,35] have been expressed in units of r . 
The interpretation of our results is not very clear, also due to the large errors. 
In any case, we want to stress that our determination can be seen at least as an 
upper limit for m ++ and m 2 ++. We expect that at the values of t/a at which we 
extracted the masses, the effects of positivity violations have already disappeared; 
this assumption is justified by our study based on perturbation theory and our 
estimation of t m i n eq. (A. 9). 

For this reason we believe that possible systematic uncertainties on the glueball 
masses could only be due to the presence of excited states and hence could affect 
our measurement only in such a way that the real values of m ++ , m 2 ++ are lower 
with respect to our determination. 

At lattice spacings a ~ 0.15 fm we notice a improvement of the RG actions with 
respect to the Wilson action; comparing with the continuum limit we find no 
significant discrepancy both for DBW2 and Iwasaki action, while for the Wilson 
action one finds 30 — 40% deviation. 

At lattice spacing a ~ 0.1 fm we find on the other hand large lattice artefacts for 
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0.5 



Iwasaki, /? = 2.5206 







2 4 6 8 

t/a 

Figure 10: The effective masses for the A± + channel, evaluated with Iwasaki action, at different 
lattice spacings. The hlled and empty squares corresponds to respectively to cq. (7.10) and 
cq. (7.12) with to = 0. For [3 = 2.5206, the filled triangles correspond to eq. (7.10) with to = a. 
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Figure 11: The effective masses for the A± + channel, evaluated with DBW2 action, at different 
lattice spacings. The hlled and empty squares corresponds to respectively to cq. (7.10) and 
eq. (7.12) with i = 0. For (3 = 0.9636, the filled triangles correspond to eq. (7.10) with t = a. 



RG actions: the result obtained with the Iwasaki action is compatible with the 
one calculated through the plaquette action at the same lattice spacing, while for 
the DBW2 action it is even further away from the continuum limit. 
If one considers our measurement as upper limit, one could conclude that the RG 
improved actions are not able to cure the problem of large lattice artefacts for the 
++ glueball mass. At very small lattice spacings one expects that the dominant 
lattice artefacts are of order a 2 ; for the Wilson action this is indeed well confirmed 
by the numerical results in the range a < 0.1 7fm, as one can see in Fig. 13. 
For alternative actions there is no reason a priori to observe the same behavior: 
while at small lattice spacings the 0(a 2 ) should in any case be the dominant one, 
at larger a it is possible that lattice artefacts are governed by higher orders (0(a 4 ) 
and higher). 

We have indeed already observed deviations from 0(a 2 ) behavior in RG actions 
for the quantity T c r and even already for the force computed at tree level in the 
previous sections. 

In Fig. 14 we report our results for r m2++; for this particular observable the 
calculation performed with the Wilson action do not show significant lattice arte- 
facts. 

At our smallest lattice spacing we do not observe a deviation from the results 
obtained with the Wilson action; one has however to notice that our errors are 
too large to make any conclusive statement. 

For our largest lattice spacing the results with the Wilson action are not available 
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Figure 12: The effective masses for the T 2 ++ channel, evaluated with Iwasaki and DBW2 
action, at different lattice spacings. The filled and empty squares corresponds to respectively to 
cq. (7.10) and cq. (7.12) with t = 0. 



22 



1 1 u 

Lollab. 


r m ++ 


r m 2 ++ 


M & P [37] 


4.21(H)(4) 


5.85(2)(6) 


GF11 [34] 


4.33(10) 


6.04(18) 


Teper [19] 


4.35(11) 


6.18(21) 


UKQCD [35] 


4.05(16) 


5.84(18) 


FP [12] 


4.12(21) 


[5.96(24)] 



Table 7: Continuum extrapolations of the two lowest glueball masses in units of r . For the 
FP action, the 2 ++ value is not extrapolated to the continuum but denotes the mass obtained 
at a lattice spacing a = 0.10 fm. 



Iwasaki action 



p 


am A ++ 


am T ++ 


2.2423 
2.2879 
2.5206 


1.11(10) 
1.20(7) 
0.72(6) 


1.50(32) 
1.17(10) 


DBW2 action 


P 


am A ++ 


am T ++ 


0.8243 
0.9636 


1.27(10) 
0.62(7) 


1.70(39) 
1.22(15) 



Table 8: Glueball masses in lattice units, Iwasaki and DBW2 action. 

(there are only know results with anisotropic lattices [37]) and it is indeed diffi- 
cult to have a reliable estimation of the mass also in our case: we decided to show 
however our results, even with the very large error bars. 

From our computation one can not deduce if for r m 2 ++ the RG actions show 
significative discretization errors and further investigations are needed to clarify 
the issue. 

In particular, the exponential error reduction proposed by Liischer and Weisz 
in [38] has been already tested for the evaluation of the ++ and 2 ++ masses 
in [39] yielding promising results. We expect that the implementation of this algo- 
rithm for RG actions could help in reducing the errors for the correlation functions 
at quite large t and hence to have more reliable estimations of the glueball masses. 
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Iwasaki action 



(3 r m A ++ (r n ) r m T ++ (r n ) r m A ++ (n) r m T ++ (n) 

2.2423 3.08(28) 3.07(28) 

2.2879 3.63(21) 4.54(97) 3.63(21) 4.54(97) 

2.5206 3.26(27) 5.31(45) 3.25(27) 5.28(45) 



Table 9: Results for r^ma for the channels + and T^ + , using the Iwasaki action. 

DBW2 action 

(3 r m A ++ (r n ) r m T ++ (r n ) r m A ++ (n) r m T ++ (n) 

0.8243 3.97(31) 5.3(1.2) 3.86(30) 5.2(1.2) 

0.9636 2.86(32) 5.62(69) 2.82(32) 5.56(68) 



Table 10: Results for romc for the channels Af + and T^ + 1 using the DBW2 action. 
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Figure 13: The ++ glucball mass normalized with ro as function of (a/ro) 2 for different 
actions. 
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Figure 14: The 2 ++ glueball mass normalized with r as function of (a/r ) 2 for different 
actions. 
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8 Discussion 



The adoption of alternative gauge actions with the purpose to reduce the lattice 
artefacts and/or to improve chiral properties in QCD simulations must be accom- 
panied with an accurate study of their properties. In particular, universality and 
scaling behavior must be tested on a large set of observables. 
We have chosen r , the deconfining temperature and the glueball masses m ++, 
m 2 ++ to perform this test for the Iwasaki and DBW2 actions, which contain the 
same operators with different coefficients. 

We computed the scale r through the force between static quarks at several values 
of (3 corresponding to the critical couplings for deconfinement and finally built the 
renormalized quantity T c r . By approaching the continuum limit we verified the 
universality between Wilson and RG action, which had been doubted in previous 
calculations where the string tension a was used to set the scale instead of r . 
Our new computation confirms that r is a more appropriate quantity, since the 
systematic uncertainties are reduced with respect to the string tension. 
From the point of view of the scaling behavior, for this particular quantity the 
Iwasaki action is able to reduce the lattice artefacts. For the DBW2 action, on 

the other hand, we observe large lattice artefacts at a ~ While they can 

be made small in r T c by choosing a r\ in the definition of the force,this choice 
leeds to large a-effects in a q g. The coefficient c\ appears to be too far away from 
its tree-level value. By investigating the force at tree level we pointed out that 
RG actions are "over-corrected" , in the sense that they can introduce even larger 
lattice artefacts than the Wilson action. This suggests that the improvement may 
work only in specific coupling ranges; moreover, it can happen that the 0(a 2 ) 
scaling violation is dominant only at very small lattice spacings, and this feature 
requests a particular care in performing the continuum extrapolation. 
Another important fact that we stressed is the violation of physical positivity 
which occur in RG actions, and more generally in actions that contain other 
terms apart from the plaquette. We clearly observed the presence of unphysical 
states in the evaluation of the static quark potential from Wilson loop correlation 
functions. 

The main problem connected with this fact is the failure in the mathematical 
assumptions which allow to apply the variational method. The method can be 
safely applied only if one manages to get rid of the unphysical modes by going at 
sufficiently large time separations - but here the statistical errors increase dramat- 
ically - or by finding a reliable procedure to project them out from the correlation 
matrix. 

We followed [27] and estimated in perturbation theory at which time separation 
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these lattice artefacts are expected to disappear and this indicative study is con- 
firmed by our numerical evaluation of the static potential. 

We stress that we did not perform a systematic study on the efficiency of the 
smearing procedure for RG actions and we applied the same criteria adopted for 
the Wilson action at analogous lattice spacings. The presence of states with neg- 
ative norm can spoil the search for an "optimal" smearing since it can mask the 
presence of excited states at small time separations, and a reasonable criterium 
must be found. 

In a second part of the work we performed numerical simulations to evaluate m ++ 
and m 2 ++ at several lattice spacings for RG actions. The main motivation was 
that the ++ glueball show large scaling deviations and hence is an ideal observa- 
tory to make a comparison with alternative actions. 

Here the violation of positivity constitutes an even stronger problem, since for 
these observables the signal is lost in noise already at small time separations. 
Due to these difficulties, our estimations of glueball masses are affected by large 
errors and hence not suitable to draw a definitive conclusion about the lattice 
artefacts; a possible solution would be to apply specific algorithms for the vari- 
ance reduction in order to extend the precision on the effective masses to larger 
values of t, where the unphysical states are supposed to be absent. 
The fact that for the Wilson action the lattice artefacts for m ++ are large and 
the glueball mass in lattice units becomes small is usually interpreted as the 
"influence" of the endpoint of the first order phase transition which for SU(3) pla- 
quette action is located in the fundamental- adjoint coupling plane at (/3/,/3 a ) = 
(4.00(7), 2.06(8)) [40]. 

This issue turns out to be very important for the next simulation with dynami- 
cal fermions; in particular, in 0(a) improved Wilson fermions it is believed that 
the clover term produces an adjoint term in the gauge action with positive cou- 
pling which could be responsible for the observed first order phase transition at 
zero-temperature in the Nf = 3 simulations [41]. The JLQCD investigations indi- 
cate that this transition is a lattice artefact restricted to strong coupling regions 
(f3 < 5.0) and it disappears if one adopts other gauge actions, like Iwasaki or 
tadpole-improved Symanzik action, or alternatively if one uses the unimproved 
Wilson fermionic action with the plaquette gauge action. 

By considering our estimation as an upper limit for the glueball masses, it seems 
that RG actions do not lead to a sensible reduction of the lattice artefacts for 
m ++, and hence the scenario remains unclear. A possible alternative to avoid 
the problem would be to adopt the usual plaquette action with a negative adjoint 
term, and studies on the scaling behavior for this kind of actions are ongoing. 
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A Positivity violation in improved gauge actions 

In [27] it has been pointed out that for lattice gauge actions containing general 
classes of loops it is still possible to construct a transfer matrix T. In particular, 
restricting to actions of kind eq. (2.1) with plaquettes and planar (1 x 2) loops, 
the transfer matrix can be defined as transition amplitude between fields config- 
urations defined on two consecutive pairs of equal time hyperplanes. But, unlike 
the Wilson case, one finds that T is no longer hermitean, and its spectral values 
A occur in pairs of complex conjugated numbers. 

Moreover, one can show that the spectrum cr(T) contains a real, positive, non- 
degenerate eigenvalue Ao such that 

A >|A|, VAeo-(T), A^A . (A.l) 

The corresponding eigenfunction may be interpreted as ground state wave func- 
tion. 

The violation of physical positivity forbids the definition of an Hamiltonian oper- 
ator; from the point of view of the spectral decomposition of connected two-point 
functions, one still observes expontential decay at large time-separations, but con- 
tributions with negative weight will appear, as remnant of the positivity violation. 
Nevertheless, the violation of positivity is a lattice artefact and one expects that 
physical positivity is recovered in the continuum limit; one should be able to de- 
fine a subspace of states with small energy with respect to the cut-off where the 
non-hermiticity of the trasfer matrix can be eliminated. 

This corresponds to the fact that for large enough time-separations, the negative 
contributions to the spectral decomposition of two-point functions should disap- 
pear. 

Following [27], one can argue that there exists a < e < 1 such that independently 
of the cutoff the following properties hold: 

(i) all spectral values A G cr(T) with |A| > eA are real and positive; 

(ii) one can define a physical Hilbert space TL p h y s G 7~t and a new scalar product 
such that O, V) new > 0, Vtf G H phys , * ^ 0. 

The existence of such e has not been rigorously proved, but it is supported for 
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example by perturbative calculations. In the next section we will estimate e by 
determining the location of unphysical poles in the free propagator. 



A.l Unphysical poles in the propagator 



In [9] the propagator associated to the action eq. (2.1) is evaluated to lowest 
order in g with covariant gauge fixing; it is defined by 

{^(x)Ai(y)) = S ij I e tk ( x - y) e l{k »- k » )/2 D^(k), (A.2) 

where A^(x) is related to the link variables by 

U{x,n) = e~ A ^ x) 
and the momenta integration corresponds to 



/ 11/ 



n/ 



We are now referring to the case L = oo, so that the momenta take continuous 

values in the Brillouin zone. 

Eq. (A.2) can be rewritten in the form 



D VT (k) = (k ) 



2\-2 



cxk v k T -\- ^ (k<j8 TU k T b av ) A TO -(/c) k a 



(A.3) 



where a is the gauge parameter, k^ = 2sin(A^/2) and is independent of a 
and has the general form 

V = % (A.4) 

with 



D = ^2 K II ^ + Yl ^hl%v{% P (lvT + q^qvp)- (A.5) 

P vj^p p>U,p>T,{p,T}U{p,v}=® 

The explicit form of f(k) can be found in [9] and will not reported here. The 
denominator of eq. (A.3) is then given by 



2\2 



A = D(k ) 



(A.6) 



(k 2 ) 2 j ^ K II q ^ + J2 ^l%M^r + Qprqup) ] , 

p v^p p>u,p>T,{p,r}U{p,u}=d 
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where q^ v for our specific form of the action eq. (2.1) takes the form 

q^ = (l-5^)[l- Cl (kl + kl)}. (A.7) 

Then eq. (A. 6) becomes 



In order to search for the poles of the propagator, one substitutes 

k = (ki, k 2 , kz, iw) 

and looks for solutions of the equation A = scanning the whole Brillouin zone. 
In general for c\ ^ one expects complex conjugated solutions 

w = *Rz(w) ± Om(«)) = w(ki, k 2 , ^3). 

Numerical investigations showed that the condition 3m(w) = (for all solutions 
at a given momentum) defines the equation f(ki, k\, k^) = which represents a 
compact 3-dimensional object. In Fig. 15 the 2-dimensional intersection with the 
plane k% = is plotted for several values of c\ corresponding to different actions. 
For our numerical studies, we discretized the Brillouin zone in finite intervals that 
could be made arbitrarily small. 

In the region enclosed by the curve one has 3m(w) = 0, while outside 3m(w) 7^ 

0. 

The Fig. 16 represents the distribution of the poles; by scanning the Brillouin 
zone and varying the intervals of the momenta, one reaches the conclusion that 
there exists a maximal value of *Rz(w), which is related to c\ by 

{ 2.063 Symanzik, tree level 
1.162 Iwasaki 
0.588 DBW2 



such that for w < w max the imaginary part 3m(w) vanishes. 

The value of w max yields an estimate of e defined above through the simple relation 

0.127 Symanzik, tree level 
i "' < 0.313 Iwasaki 
0.555 DBW2 
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Figure 15: The curves defined by the condition 3m(io) = in the ^3 = plane. Inside the 
curve 3m(w) = 0, while outside 3m(w) ^= 0. 
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Figure 16: Distribution of the poles for different actions, scanning the Brillouin zone. 
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Evaluating observables from exponential decays of correlation functions, one ex- 
pects the presence of unphysical states unless t ^> t m i n , where t m i n in lattice units 
is given by 

{ 0.484 Symanzik, tree level 
0.860 Iwasaki (A.9) 
1.702 DBW2 

One has to remember that our estimation of the unphysical poles has been per- 
formed in perturbation theory and hence has to be considered as an indicative 
evaluation. 

B Numerical results 

The tables 11 and 12 report the data for the potential and the force evaluated 
at finite lattice spacing with Iwasaki and DBW2 actions. r\ defines the tree- 
level improved force and is defined in eq. (4.7). For the potential only the naive 
definition is reported. 
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Table 11: The potential and the force in lattice units for the Iwasaki action. 





ri/a 


a 2 F(rj) 


r/a 


aV(r) 


2.1551 


1.4858 


0.4059(21) 


2 


0.9519(20) 




2.4626 


0.2965(75) 


3 


1.2487(61) 




3.5663 


0.2585(98) 


4 


1.515(12) 




4.6059 


0.249(29) 


5 


1.783(12) 


2.2879 


1.4858 


0.6055(29) 


2 


0.7730(12 




2.4626 


0.19796(76) 


3 


0.9707(33) 




3.5663 


0.1689(70) 


4 


1.1418(41) 




4.6059 


0.1622(64) 


5 


1.3047(95) 




5.6082 


0.155(10) 


6 


1.460(16) 




6.6036 


0.148(16) 


7 


1.608(33) 


2.5208 


1.4858 


0.21301(25) 


2 


0.60342(29) 




2.4626 


0.11734(39) 


3 


0.72099(46) 




3.5663 


0.09011(70) 


4 


0.8107(14) 




4.6059 


0.0803(10) 


5 


0.8902(20) 




5.6082 


0.07575(61) 


6 


0.9658(36) 




6.6036 


0.0744(32) 


7 


1.0378(37) 




7.5977 


0.0711(15) 


8 


1.1098(74) 


2.7124 


1.4858 


0.17512(12) 


2 


0.52165(14) 




2.4626 


0.08691(14) 


3 


0.60856(20) 




3.5663 


0.06138(53) 


4 


0.66994(56) 




4.6059 


0.05194(81) 


5 


0.7219(13) 




5.6082 


0.04693(85) 


6 


0.7688(20) 




6.6036 


0.0451(23) 


7 


0.8139(11) 




7.5977 


0.0425(10) 


8 


0.8564(10) 




8.5915 


0.0420(15) 


9 


0.8984(23) 




9.5854 


0.0414(19) 


10 


0.9398(25) 
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Table 12: The potential and the force in lattice units for the DBW2 action. 



a 
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77 


T?( r \ 
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r / a 


a v \r ) 
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1 1 79KfQQ N \ 
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1 /ii n/1 1 "\ 
1.41U(11 J 
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n 9/19/1 fCQ"\ 
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1.00o(ZU ) 




0. ( OOO 


U.Z0U(ZU J 
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i .oy4(oo j 


U.oz4o 


1 OOQQ 

l.yzoo 


U.zyoy i (Ul ) 


o 
z 


U. f 00 1 (11) 




O 77QQ 


n 1 QQ7f 1 9^ 
U. loo i (1^) 


Q 

o 
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0. ( 400 


U. 1U10(1Z ) 


/I 
4 


l.UoO 1 (1 I ) 




/I 7QK/I 
4. / oU4 


U. looz(Uz ) 





l./ooy( I U) 




0. i 000 


U. 14 I K){i) 1 J 





1.0oD(10 ) 




6.7843 


0.1494(61) 


7 


1.536(2U) 


u.yooo 
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U.zUbTU(z ( ) 


z 


U.oooy3(o4 ) 




Q 77QQ 


U. 110U(1U ) 
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o 7/1 qc 
0. ( 400 
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/I 
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/I 7QK/1 
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U.U / UU(o ( ) 


7 


U.y / < (1U) 




( .olZU 


U.U04U^4 I J 


O 

o 
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1.04 


1.9233 


0.18182(26) 


2 


0.49710(24) 




2.7793 


0.09392(54) 


3 


0.59094(58) 




3.7435 


0.06847(84) 


4 


0.65914(76) 




4.7364 


0.0575(18) 


5 


0.7180(26) 




5.7535 


0.0544(12) 


6 


0.7724(21) 




6.7843 


0.0527(22) 


7 


0.8254(42) 




7.8120 


0.0498(23) 


8 


0.8753(65) 




8.8296 


0.0481(42) 


9 


0.9230(51) 




9.8369 


0.0485(44) 


10 


0.9713(56) 




10.8365 


0.0464(34) 


11 


1.0175(95) 




11.8310 


0.048(13) 


12 


1.067(20) 
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